Self-consistent study of electron confinement to metallic thin films on solid surfaces 
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We present a method for density-functional modeling of metallic overlayers grown on metallic 
supports. It offers a tool to study nanostructures and combines the power of self-consistent pseu- 
dopotential calculations with the simplicity of a one-dimensional approach. The model is applied 
to Pb layers grown on the Cu(lll) surface. More specifically, Pb is modeled as stabilized jellium 
and the Cu(lll) substrate is represented by a one-dimensional pseudopotential that reproduces ex- 
perimental positions of both the Cu Fermi level and the energy gap of the band structure projected 
along the (111) direction. The model is used to study the quantum well states in the Pb overlayer. 
Their analysis gives the strength of the electron confinement barriers at the interface and at the 
surface facing the vacuum. Our results and analysis support the interpretation of the quantum well 
state spectra measured by the scanning tunneling spectroscopy. 



I. INTRODUCTION 

The ability to manufacture desired nanostructures on 
solid surfaces by controlled growth or atomic manipu- 
lation techniques has increased enormously during the 
last decadesi^ At the same time, spectroscopic meth- 
ods to study different physical properties and phenomena 
of these structures have also experienced a huge devel- 
opment. In understanding and interpreting the ensuing 
experimental results rich of quantum phenomena, accu- 
rate theoretical and computational modeling plays a vital 
role. 

A vifidely studied phenomenon is the growth of thin 
Pb films or extended Pb islands on solids, for exam- 
ple, on Si or Cu surfaces. The growth of Pb provides 
a laboratory to test the the so-called quantum size ef- 
fects (QSE) arising due to the electron confinement per- 
pendicular to the surface. The confinement results in 
discrete energy levels, the so-called quantum well states 
(QWS's). With the increasing film thickness new QWS's 
become occupied producing oscillations in the total en- 
ergy, work function and other physical properties. The 
oscillations in the energy are the origin of the "magic" 

islands heightS|3.4,5i6i7,8,9,10,ll,12,13,14,15,16 

The most important feature characterizing an over- 
layer or a nanoisland is its height. However, measure- 
ments based on different physical processes may provide 
different values. The x-ray diffractionii can be used to di- 
rectly determine the number of atomic monolayers (ML) 
in thin films. The scanning tunneling microscope (STM)^ 
gives the thickness of finite nanoislands. In the helium 
atom scattering (HAS) the specular returning point of He 
atoms is measured."*'^ - STM and HAS reflect the electron 
density profile of the surface.^' The scanning tunneling 
spectroscopy (STS) measures the QWS spectrum and the 
width of the quantum well confining the electrons is eval- 
uated from it.^^ In addition to the overlayer- vacuum sur- 
face profile, the QWS spectrum is sensitive also to prop- 
erties of the substrate-overlayer interface. The informa- 



tion about the confining quantum well is important not 
only for the determination of the the film thickness but 
also for the understanding and controlling the so-called 
" electronic growth mode" . 

In this work we model thin Pb films growing in the 
[111] direction on the Cu(lll) surface. We focus is on de- 
termining the surface and interface barriers of the confin- 
ing quantum well and the ensuing QWS's. Pb/Cu(lll) 
was recently studied by Otero et. aliA& They fitted 
the QWS spectra measured by STS by using the finite- 
potential-well model and unexpectedly found that actu- 
ally the infinitely high potential barriers give a better 
fit than the more realistic finite barriers. One of our 
main aims is to resolve this dilemma. There are also re- 
cent STS measurements of QWS's in the Pb layers on 
Si(lll)»ii*SL2i In these studies the overlayers are thin- 
ner and the QWS energy range is smaller than in the 
Pb/Cu(lll) measurements by Otero et. a,lM> The effec- 
tive thickness of the Pb overlayer is then obtained by 
fitting the energy differences of the consecutive states. 
This is not very accurate and there are contradictory re- 
sults for the thickness of the wetting layer, ranging from 
1 ML7'8'9 to 3 ML.i"."i20,2i 

There are several theoretical works based on 

the density functional theory (DFT) devoted to 

study QSE in thin films. First-principles atom- 
istic approachea^^i^3,^4i^5,26,27,28,29^30 or jeiii^jj^ 

j^Qjjgjg3ij32j33j34j35 Y^g^yQ been used. Usually, due to 
the high computational demand the substrate is not 
included, but the electronic structure is calculated for 
a slab describing the overlayer. The work by Hong et 
al.'^^ for Pb/Si(lll) is one exception. On the other 
hand, there are several simple analytic models for the 
confinement barriers ^2i22i^ but a priori assumptions 
about the barrier type, its position or quantum numbers 
of the measured states can lead to an erroneous inter- 
pretation of the experiments. The simplifications made 
in the modeling hinder effectively the analysis of the 
experimental results. 
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In the present work we perform self-consistent elec- 
tronic structure calculations by modeling the Cu(lll) 
substrate with a one-dimensional (ID) pseudopotential 
and the Pb overlayer by stabilized jcUium. The model 
includes the effect of the substrate-film boundary so that 
the penetration of the QWS into the substrate is realis- 
tically described. Our self-consistent results allow us to 
also study the above-mentioned simple analytical mod- 
els and point out their deficiencies as well as the most 
important factors for the proper description of the elec- 
tron confinement. This knowledge is specially important 
when using these models in analysing the STS results for 
completely covered substrates or for systems with wet- 
ting layers of unknown thicknesses. 

In a previous publication^^ we applied successfully 
the ID-pseudopotential - stabilized-jellium model to gain 
physical insight into the "magic" heights of Pb islands 
on Cu(lll) up to 23 ML of Pb,^ A pragmatic aim of 
the present paper is to document the construction of 
the unscreened ID-pseudopotential and provide a sim- 
ple parametrization which can be used in future studies, 
e.g. for different nanostructures on surfaces. 

The rest of the paper is organized as follows: in Sec. 
nil we report the practical steps for the construction of 
the Cu(lll) pseudopotential. Then, we describe ana- 
lytical models for the confinement barriers, whose reli- 
ability is analyzed by applying them to the results of 
the self-consistent calculations. Section Hill demonstrates 
the resulting self-consistent electronic structures for the 
Pb/Cu(lll) system. The confinement barriers are ana- 
lyzed and an improvement to the theory is introduced in 
order to acquire accuracy in their determination. Sec- 
tion IIVI contains the conclusions of the work. Atomic 
units {i.e., h — e'^ — m — 1 and distances measured in 
Bohr radius units ao = 0.53 A) will be used throughout 
this work, unless otherwise specified. 



II. THEORY 
A. The model 

The present calculations are performed in the frame- 
work of the DFT— within the local density approxima- 
tion (LDA)^2i^ and jellium-type models. Instead of finite 
Pb islands we consider infinitely-extended films on the 
Cu surface. This is justified because in the experiments 
considered the character istic l ateral dimension of the Pb 
islands is around lOOOao '3*191 so that the lateral electron 
confinement effecs are vanishingly small. Moreover, we 
assume the perfect translational invariance, i.e. a homo- 
geneous free electron gas, along the surface {xy plane). 
We use a jellium model for the Pb overlayer and for the 
Cu substrate we construct a pseudopotential which varies 
only along the z direction perpendicular to the surface. 
Then the Kohn-Sham equations have to be solved numer- 
ically only in the z direction (ID problem). This enables 
the calculation of electron wavefunctions extending deep 



into the Cu substrate and the modeling of systems having 
tens of Pb ML's. 

The Kohn-Sham equations are discretized in a regular 
one-dimensional point mesh and solved with the Raylcigh 
quotient multigrid method^^*^ implemented in the real- 
space MIKA package^ for electronic structure calcula- 
tions. Hence, single-particle wave-functions are taken to 
be of the form 

<f(r) = ^„(z)e*ll--|l, (1) 

where ipniz) is the wavefunction in the direction perpen- 
dicular to the surface, and plane waves are used for the 
surface parallel directions. The eigenenergies are given 
by 

II 

enfc|| = + y, (2) 

where e„ is the eigenvalue of the n-th perpendicular state 
ipn{z). The eigenvalue e„, obtained self-consistently, is 
the bottom of the n-th subband. For finite and periodic 
systems in the z direction Dirichlet and periodic bound- 
ary conditions are used, respectively. 

The effective or screened potential of the Kohn-Sham2& 
equations in the z direction is written as 

VM= [ "-(^0 - ^ y.e[n-(z)] + Vp,{z), 

J \z-z\ 

(3) 

where the first term on the right-hand side is the Hartree 
term Vh{z), which includes the electron density n_(z) 
and the neutralizing rigid positive charge density n+{z). 
The second term gives the LDA exchange-correlation po- 
tential. The third term accounts for the pseudopotential 
that improves the simple jellium scheme. For the sup- 
ported overlayer system Vps{z) has two contributions, one 
from the Pb overlayer and the other from the Cu(lll) 
substrate. The free electron-like character of Pb at the 
Fermi level justifies the use of the stabilized jellium or 
averaged pseudopotential^^i^ approach to model Pb. In 
practice, the jellium model allows us to simulate any Pb 
overlayer thicknessi^i The Pb contribution to Vps{z) sta- 
bilizing the electron gas at the density corresponding to 
Ts = 2.30ao is a constant shift Vstab relative to the vac- 
uum level and restricted in the region of the positive 
background charge. The stabilized Pb provides a proper 
work function so that the spilling of the electron density 
into the vacuum is well described. It also gives a proper 
value for the bottom of the valence electron band. This 
guarantees the correct Fermi wave length \p, which is 
of crucial importance for the properties related to the 
electron confinement. 

For the Cu(lll) substrate it is necessary to use a pseu- 
dopotential which accounts correctly for the width of the 
energy gap and its position with respect to the Fermi 
level. In this way, we obtain the correct confinement po- 
tential also at the Cu(lll)-Pb interface. To obtain the 
ID-pseudopotential we use the procedure described in 
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the following subsection. We start from a model poten- 
tial and build an unscreened pseudopotential for Cu(lll) 
in two steps. 

B. Generation of the Cu(lll) ID-pseudopotential 

1. Bulk pseudopotential (V^ step) 

Chulkov et al^ proposed a fully screened ID-model- 
potential which varies only in the z direction perpendic- 
ular to the surface. The model potential is successfully 
used to study, for example, the dielectric response func- 
tion and lifetimes of excited states4SiiLi& The crucial 
point here is the proper description of the energy band 
gap and work function. Moreover, the wave functions 
are correctly described not only outside the substrate, 
but also inside it. This is an important ingredient in the 
present application. 

Because in this work we cover the Cu(lll) surface 
with several ML's of Pb we skip the surface part of the 
ID-pseudopotential. Nevertheless, it is also possible to 
build a pseudopotential which reproduces the surface and 
image states4^ The bulk oscillating function of the ID- 
model-potentialiS, is 

f 2tt \ —d d 

Vmodel[z) = AlO+Ai COS I —Z 1 , — < Z < - [A) 

where d = 3.943ao is the interlayer spacing in Cu in 
the [111] direction and Aiq and Ai are fitting param- 
eters. Using periodic boundary conditions at ±d/2 the 
Kohn-Sham equations are solved for the fixed Veff{z) — 
Vmodei{z) potential. With the eigenfunctions obtained 
and with the experimental work function we calculate 
the electron density profile. Integrating over z we obtain 
the mean density with rg — 2.55ao. This value is close to 
the experimental = 2.67ao for 4s electrons. 

Once we have computed the density, it is straightfor- 
ward to obtain the corresponding Vxc{z) potential. It 
is more challenging to calculate the Hartree Vh{z) term 
because in absence of vacuum we don't know the height 
of the dipole barrier, and therefore, where to fix the en- 
ergy origin for the Hartree term inside the bulk. To solve 
the problem we fix provisionally the zero of the Hartree 
potential at the boundaries of the periodic cell. After 
adding a homogeneous neutralizing positive background 
of Ts = 2.55aoi the Hartree potential can be evaluated. 
Finally, we can calculate from Eq. (O the unscreened 
and periodic pseudopotential Vps{z) by subtracting from 
the effective potential the Vh{z) and Vxc{z) terms. 

We have obtained, in this first step, the unscreened 
pseudopotential for periodic bulk calculations. Perform- 
ing a self-consistent calculation for the bulk Cu(lll) with 
this unscreened pseudopotential Vps{z) and a positive 
background density corresponding to = 2.55ao, we 
recover the experimental Fermi level with respect to the 
vacuum. 



2. Pseudopotential for a semi-infinite system (2""^ step) 

The unscreened ID-pseudopotential from the first step 
is suitable for bulk calculations. But it cannot be used in 
slab calculations because the zero of the Vh{z) potential 
was arbitrarily chosen in the previous step. We complete 
the pseudopotential for slab calculations in this second 
step. 

We build a semi-infinite Cu(lll) slab by repeating the 
pseudopotential for one ML of Cu. The slab is thick 
enough to avoid interaction between surfaces and finite- 
size effects in determining the band structure. Then, 
enough vacuum is added on both sides to annul bound- 
ary effects at the borders of the calculation volume. In 
the self-consistent calculation for this system the electron 
density spills out of Cu(lll) to the vacuum giving rise 
to the dipole Coulomb barrier which places the Fermi 
level (and the whole band structure including the energy 
gap) at the incorrect position with respect to the vacuum 
level. 

To correct the work function we shift the potential by a 
constant inside the Cu slab. Here we define the Cu(lll) 
edge to be one half ML above the last atom plane. But 
the potential shift also changes the electron spilling into 
the vacuum and the dipole barrier. Thus, we have to find 
the potential shift iteratively so that the correct value 
for the work function is recovered. The value of the shift 
depends also on the possible smoothing (see below) of 
the positive background charge profile at the surface. 

The ID-pseudopotential vanishes suddenly at the 
Cu(lll) surface, producing a discontinuity in the total 
potential. We have examined different ways to obtain a 
pseudopotential without the discontinuity. However, the 
final results obtained do not differ from the results ob- 
tained with the scheme described above and presented in 
this paper. 

In summary, we have obtained a ID-pseudopotential 
which with the positive background density of = 
2.55ao reproduces correctly the experimental Cu(lll) 
work function and the [lll]-projected band structure in- 
cluding the band gap. The method presented here for 
the Cu(lll) pseudopotential generation is extensible to 
other substrates as well. 



3. Analytic expression of the ID-pseudopotential 

The unscreened pseudopotential obtained can be fitted 
using the form of Eq. 2]for the screeneed model-potential. 
The result is, = -1.89 eV and Ai = 5.01 eV . 

This pseudopotential is obtained with the positive charge 
background density corresponding to — 2.55ao and a 
surface profile of the positive charge given by a Fermi-like 
distribution 

= Jo: ■ (5) 
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Above, Dq is the surface edge position and Az = 0.09ao 
accounts for the smoothing at the edge. The smoothing 
is used for numerical reasons. 



Phase accumulation model for confinement 
barriers 



We use the onedimesional-pseudopotential described 
above to simulate the Pb/Cu(lll) system. In order to 
rationalize the results in a quantitative way, we employ 
simple analytical models reducing the information about 
the confinement barriers in few parameters. In particu- 
lar, we are going to use the so-called phase accumulation 
model. ■^^ In this model a QWS in a confining potential 
fuUfils the equation 



2KD 



l^Cu(lll)-Pb 



^Pb~vac 



2vr(n~l), (6) 



where kz is the wavevector corresponding to the QWS ki- 
netic energy, D is the width of the potential well to model 
the Pb film, (ppb-vac and 0cm(iii)-P6 a-re the phases of 
the eigenfunction accumulated at the vacuum surface and 
at the interface, respectively, and n — 1,2,3,... is the 
quantum number of the QWS. For an infinitely deep 
square potential well, the phase accumulated on each 
surface is — tt. For softer potential barriers the value in- 
creases. 

The phase shifts, usually calculated for model poten- 
tials given analytically, incorporate the physical features 
of the QWS's in a simple manner. An equivalent but 
more intuitive way is to apply the infinitely deep square 
potential well so that (jipb-vac = (j^Pb-Cu = — tt, but us- 
ing an effective width of D' = D + 5cu{iii)~Pb + ^Pb-vac- 
Then Eq. 10 gives 



<^Cti(lll)-Pb + 5pb-va 



Tm 



D. 



(7) 



Here, Scu{iii)-Pb and 5pb-vac arise from the wavefunc- 
tion penetration into the Cu(lll) substrate and the spill 
out into the vacuum, respectively. The idea is that D' 
should give a reliable estimate of the actual thickness of 
the overlayer. Thus, as the phase shifts, also the effective 
well width depends on the QWS eigenenergy £„. It has 
been shown that the energy spectrum is very sensitive to 
the positions of the barriers and relatively insensitive to 
the barrier height ji^ Therefore, a mean surface shift 5q 
can reproduce quite accurately the spectrum. According 
to Eqs. © and Q the surface position shifts and phase 
shifts are related by 



(8) 



Iky, 



We consider below two analytical expressions for the 
surface- vacuum phase shift, the finite-potential-step and 
the image-potential models. The former gives the energy- 
dependent phase shif t^^i^" 



^ Pb—vac 



— 2 arctan — 



h 



(9) 
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FIG. 1: Electron density (solid line) and average positive 
background density (dashed line) for 4 ML's of Pb on the 
Cu(lll). The origin « = is at the Cu-Pb interface. The 
dotted curves show the electron densities of the free-standing 
Cu(lll) and Pb slabs. The dash-dotted curve shows the 
charge transfer in the combined system relative to the free- 
standing Cu(lll) and Pb slabs (Ap = pp^/cuiiii) ~ Pcu(iii) ~ 
pp^). The inset sketchs the symmetric geometry used in the 
calculations. 



where kyac = •\/2(— e) with the energy eigenvalue e mea- 
sured with respect to the vacuum level. This model has 
been used for the Cu-Pb interface barrier, tooii^ The 
image-potential model gives the phase shifii^ 



i'Pb-vac — TT 



3.4 eV 



(10) 



The phase shift corresponding to the QWS wavefunc- 
tion penetrating into the Cu(lll) depends on the posi- 
tion of the QWS energy eigenvalue relative to the energy 
band gap. For example, the empiric formula^i 



4>Cu(iii)-Pb = 2arcsin 



e„ - e. 



(11) 



has been used. Above, and are the upper and lower 
edges of the band gap, respectively. 



III. RESULTS AND DISCUSSION 
A. Electronic structure 

The main results of this paper are obtained with the 
pseudopotential generated in the previous section for the 
Cu(lll) surface and the stabilized jellium model for the 
Pb overlayer. In the inset of Fig. ^ we sketch the ge- 
ometry used in the calculations. In the middle there are 
25 ML's of Cu and Pb overlayers corresponding up to 23 
ML's (1 ML of Pb is 5.41 oo) are attached symmetrically 
on both sides. Complementary calculations have been 
performed for unsupported Pb slabs, too. 



5 



(a) (b) 




FIG. 2: (a) Effective potential 14// (z) (solid line) and its components, the exchange-correlation potential Vxc{z) (long-dashed 
line), Hartree term Vh{z) (short-dashed line) and Pb stabilization potential Vstab (dotted line) for 4 ML's of Pb on the Cu(lll) 
surface. The origin z = is at the Cu-Pb interface and the energies are measured with respect to the vacuum level. The 
dash-dotted curve is the potential obtained by using the stabilized jellium model also for Cu. The straight horizontal line 
denotes the Fermi level ef- The vertical dotted lines represent the surfaces of the 4 ML Pb slab and the diamonds mark the 
position of atomic Cu planes. The shadowed regions distinguish the energy and space regions where localized QWS's can exist, 
(b) Energy eigenvalues as a function of the wavevector k^. The open circles represent states extended over the whole system 
while the filled circles correspond to QWS's trapped in the Pb layer. The dotted parabola is the free-electron-model band in 
Pb plotted as a guide to the eye. The dark gray and ligh gray areas correspond to the Cu(lll) energy gap and the potential 
well in Pb, respectively. The Cu and Pb free-electron-model Fermi wavevectors are marked for reference. 



Fig. ^shows the electron and positive background den- 
sities corresponding to the coverage of 4 ML's of Pb. We 
notice in the electron density in Pb six Friedel oscillations 
with the wavelength of about half Fermi wavelenght. In 
the present case they are produced by the 6th QWS band, 
the bottom of which has dropped well below the Fermi 
level. The 7th QWS band is still unoccupied. This kind 
of well-developed Friedel oscillation pattern commensu- 
rate with the thickness of the overlayer is characteristic 
for the stable "closed shell" or "magic" overlayer sys- 
tems. In Cu(lll) the electron density oscillates strongly, 
according to Fig. ^ as a consequence of the oscillating 
pseudopotential. 

Fig. n gives also the electron densities of the corre- 
sponding free-standing Cu(lll) and Pb slabs as dotted 
curves. It is remarkable that the electron density moves 
from Pb towards Cu(lll), indicating a possible increase 
in the effective width of the Pb slab. The charge trans- 
fer obtained when comparing the Pb/Cu(lll) system 
with its free-standing counterparts (Ap — Pp^/cu{iii) ~ 
/^cu(iii) ~ Ppb) reveals charge accumulation at the inter- 
face. A small amount of charge (~ 5% of one Pb ML 
charge) is transferred from Pb to Cu in order to equili- 
brate their chemical potentials, ~ —4.1 eV and —4.94 eV 
for Pb and Cu, respectively. The charge transfer oscil- 
lates slightly as a function of the Pb slab thickness due 
to the oscillations in the Pb/Cu(lll) Fermi level. These 
oscillations are in turn a consequence of the electron z- 



confinement in the Pb overlayer. 

In Fig. Ufa) the effective potential Vg// and its com- 
ponents Vxc, Vh and Vstab are shown; the total pseu- 
dopotential Vps is not plotted for clearness. The effective 
potential corresponding to the same system, but modeled 
with the stabilized jellium for Cu, is also drawn (dash- 
dotted line), for comparison. Notice that it gives approx- 
imately the mean value of the Cu(lll) pseudopotential. 
The dark gray shadowing gives the rough spatial exten- 
sion where the localized QWS's can appear because of 
the Cu(lll) energy gap. The light gray area marks the 
potential well between the vacuum barrier and the Pb-Cu 
interface barrier. This well is created because the poten- 
tial in Pb is '--^ 2 eV lower than the average Cu potential 
(notice the dash-dotted curve). 

In Fig. |2Ib) the eigenvalues are plotted as a function 
of the wavevector fc^. Actually, as will be explained be- 
low, it is not a straightforward task to assign a wavevec- 
tor kz to each state. Again, the dark gray area marks 
the position of the energy gap induced by the Cu(lll) 
potential and the light gray area shows the lowest ly- 
ing QWS's bound by the Pb-Cu interface barrier. The 
QWS's appearing in these energy regions are plotted with 
filled circles (Fig. [Slshows an example of the correspond- 
ing eigenf unctions (solid line)). The QWS's fall on the 
parabola (dotted line) which is the free-electron-model 
band for Pb. The states extending over the whole Cu-Pb 
system and forming a continuous band are plotted with 
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open circles. The nearly parabolic band reflects the free- 
electron character in our description of the Cu substrate. 
The QWS's in the lower, light gray area are supposed to 
have a minor relevance on the electronic properties of the 
system. In contrast, QWS's in the upper dark gray area 
play an important role, because increasing the Pb thick- 
ness lowers the QWS energy and they become occupied 
one by one. 

The repeated occupation of new QWS's as a function 
of the overlayer thickness produces oscillations in all the 
electronic properties. In particular, oscillations in the 
energy give rise to Pb overlayers of especial stability. We 
studied recently these "magic" island heights using the 
present model for the electronic structuresii^ 

As pointed out above, we divide the eigenstates into 
QWS's and bulk states. We use different methods to 
assign wavevectors for these two groups of states. For 
bulk states, of minor importance in this study, we use 
the particle-in-a-box wavevector kz = irn/D, where n is 
an integer and D the box width. On the other hand, 
QWS's are very localized in the Pb slab and we obtain 
their wavevectors from the confinement kinetic energy 
as kz — \/2ek. is the energy difference between the 
QWS energy eigenvalue and the average potential value 
around the center of the Pb slab. The wavelengths can be 
measured also from the eigenfunctions and similar values 
are obtained. Below, the kz vector turns out to be a 
relevant parameter when evaluating the phase shifts at 
the interfaces. 



B. Determination of confinement barriers I 

1. Comparison with the experiments 

Scanning tunneling spectroscopy is capable of resolving 
energy levels of QWS's. Nevertheless, the interpretation 
of the experimental results is still difficult The results 
obtained by our model provide a good tool to enlighten 
the problems encountered and to suggest a simple but 
realistic picture with relevant parameters. 

Fig. O shows the calculated eigenvalues as a function 
of the number of Pb ML's. The comparison with the 
QWS energy levels measured by Otero et al^ shows a 
good agreement for large coverage heights. Below 6 ML's 
the correspondence is a bit worse. This can be due to 
an inaccurate description of the Cu-Pb interface or due 
to the interaction with the Cu d electrons omitted in 
the calculations. The importance of these phenomena 
decreases as the number of ML's increases. 

To explain the measured energy eigenvalues Otero et 
alm^ tried first the finite-square-potential-well model but 
the fit was not satisfactory. However, they obtained 
a very good agreement by using the infinite-square- 
potential-well model or by decreasing the width of the 
finite square well by ~ 2.6ao. The fact that the infi- 
nite potential well produces much better results than the, 
in principle, more realistic finite potential well is coun- 
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FIG. 3: Energy eigenvalues of QWS's in the Pb overlayer 
on the Cu(lll) surface as a function of the number Pb ML's. 
Open circles and filled diamonds give the calculated and ex- 
perimental STS^' values, respectively. The dotted curves link 
points calculated with the infinite potential well model (See, 
Sec. irm with the effective well width of D' = D-|-4.65ao. D 
is the nominal width given by the number of Pb ML's. The 
energies are given with respect to the Fermi level. 



terintuitive, because it has been proven for slabs^^iiS 
wires^^i^ and clusters^ that the real potential profiles 
are soft. In fact, it can be seen in Fig. that the 
self-consistent potential is soft with an effective width in- 
creasing as a function of the QWS energy level. On the 
other hand, this effect seems not to be counterbalanced 
by the confining potential on the Cu-Pb interface. The 
QWS's penetrate inside the Cu(lll) (see Fig. ]^ and the 
electron density in Fig. Amoves towards Cu. Neverthe- 
less, a very good agreement with the experiments is ob- 
tained with both the infinite-square-potential-well model 
and soft self-consistent potentials. In addition, in Fig. 
13 we show the eigenvalues (linked by dotted curves) ob- 
tained with the infinite-square-potential-well model but 
by using well widths increased by 4.65ao in order to take 
into account the wavefunction spill-out. As can be seen, 
a good fit is obtained also in this way. 



2. Calculation of the phase shifts 

In order to clarify the reasons for these contradic- 
tions or dissimilar results, we use the phase- accumulation 
model and evaluate the phase shifts for the confinement 
barriers in our self-consistent calculations. We use the 
unsupported Pb slab to calculate the phase shift 4>pb-vac 
and then we can evaluate the phase shift (t>cu{iii)-Pb in 
the Pb/Cu(lll) system. 

The procedure used is to choose a QWS with an en- 
ergy just below the vacuum level and then to identify the 
corresponding quantum number n. After the wavevec- 
tor kz is calculated from the kinetic energy e = k1/2 
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as explained at the end of Section IIII Al the phase shift 
4>Pb-vac is evaluated for the QWS energy level by Eq. 
(jHJ. Increasing the slab thickness, the energy eigenvalue 
decreases and we calculate 4>pb-vac{£) as a function of the 
energy. The phase shift depends on the potential profile 
of the surface. Therefore we check that it remains unal- 
tered for thicknesses over ~ 2 ML's. Once 4>pb-vac{£) has 
been calculated by using the unsupported Pb slab, the 
phase shift 4>cu(iii)~Pb{£) can be evaluated in the same 
way using the results of the Pb/Cu(lll) system. Simi- 
larly to the above procedure for the phase shifts, we also 
determine, according to Eq. {Tj) , the shifts in the surface 
and interface positions 5pb-vac(f) and (5cu(iii)_p6(e), re- 
spectively. 

The identification of the quantum number n is easy for 
states in free-standing Pb slabs. But it is difhcult for the 
QWS's in the Pb/Cu(lll) system because, in principle, 
it is not known how many states are hidden as resonances 
in the region of delocalized states in Cu (the white area 
between the gray ones in Fig. Olb). Nevertheless, the 
local DOS integrated over the Pb overlayer reveals the 
number of the QWS resonances and, in addition, we have 
the possibility to plot the eigenfunctions in order to check 
the identification. 

The dependences of the phase shifts on the energy 
eigenvalue are shown in Fig. 0] We notice that nearly in 
the whole energy range the phase shift at the Cu(lll)- 
Pb interface is larger than that at the vacuum side. This 
means that the wavefunction penetration into Cu(lll) is 
larger than the spill out into the vacuum. The behaviour 
of the phase shift 0pb_^ac(e) can be compared with that 
obtained in the three-dimensional DFT calculations with 
more realistic atomic pseudopotentials. For example, the 
results by Wei and Chou'^^ are in a qualitative agreement 
with the preset ones. 

The phase shifts of the image-potential model, Eq. 
()10|l . reproduce quite well our results with the exception 
of the energies close to the vacuum level. This is due 
to the exponential decay of the DFT potential into the 
vacuum. The empirical Eq. gives good results at 

the Cu-Pb interface if one corrects for the known down- 
ward shift of TTi^ The finite-square-potential-well model 
does not produce phase shifts in agreement with the self- 
consistent results. 

We plot in Fig. 01b) the shifts of the barriers in the 
infinite-square-potential-well model according to Eq. Q 
as a function of the energy eigenvalue. In the energy 
window from -1 to 3 eV accessible in the STS measure- 
ments the effective width increases more at the Cu(lll)- 
Pb interface than at the Pb-vacuum surface. The shift of 
the effective Pb-vacuum barrier is in agreement with the 
jellium-model calculations by StrattonS who obtained 
the average shift of 1.41ao between the bottom of the 
potential well (at -9.5 eV) and the Fermi level. 

Combining the results of Fig. ^and Eqs. ® or ((TJ, 
the eigenvalues in Fig. 13 are reproduced. (The energy 
eigenvalues have to be shifted downwards by the width 
9.5 eV of the occupied energy band.) A good qualitative 
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FIG. 4: (a) Cu(lll)-Pb and Pb-vacuum phase shifts as a 
function of the QWS eigenenergy given with respect to the 
Fermi level. The solid curves give the results of our DFT 
calculations. The dashed curves denote the values obtained 
by using Eqs. IIOII and IIH for the Pb-vacuum and Cu(lll)- 
Pb phase shifts, respectively (vr is added to Eq. ||TTJ)- The 
dotted curve corresponds to the finite-potential-step model for 
the Pb-vacuum barrier, (b) Shifts in the effective Cu(lll)- 
Pb interface and Pb surface positions are given by solid and 
dashed curves, respectively. 



behaviour is obtained simply by using the mean effec- 
tive width D' which is the ideal width D increased by 
~ 4.65ao. The dotted curves in Fig. El link the eigenen- 
ergies calculated in this way. 

We compare in Fig. Oour self-consistent DFT calcula- 
tion for the Pb/Cu(lll) system with two simple model 
calculations. There are 5 ML's of Pb on the Cu(lll) 
surface and we plot the electron density of the QWS at 
£ = 0.53 eV (solid line). The corresponding state in a 
free-standing Pb slab (dashed line) is obtained by shift- 
ing the left Pb-vacuum boundary by 0.75ao to the left 
in order to mimic the larger penetration of the wavefunc- 
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FIG. 5: Electron density of a QWS in the 5 ML thick Pb 
overlayer on Cu(lll). The solid curve results from the self- 
consistent DFT calculation for the Pb/Cu(lll) system. The 
dashed and dash-dotted curves are obtained using the free- 
standing Pb slab and the infinite-square-potential-well model, 
respectively. The gray region denotes the Pb overlayer on 
Cu(lll). For the details, see the text. 



tions into Cu(lll) than into vacuum. The value of 0.75ao 
reproduces the correct infinite-barrier shift in Fig. Efb). 
The eigenenergy of this Pb slab state is also 0.53 eV. 
On the right surface, of course, both wavefunctions over- 
lap and at the left surface the slab wavefunction tries to 
mimic, without oscillations, the decay of the wavefunc- 
tion in the Pb/Cu(lll) system. The corresponding state 
calculated using the infinite-square-potential-well model 
with the appropriate shifts of the barriers (1.66ao and 
2.37ao from Fig. EJb)) is also given in Fig. \E\ (dash- 
dotted line). Note that this simple model gives the cor- 
rect energy eigenvalue and describes well the wavefunc- 
tion inside Pb but it cannot account for the decay into 
the Cu substrate. In conclusion, the two simple models 
reproduce quite well the eigenfunctions inside Pb and at 
the Pb- vacuum barrier, but not beyond the Cu(lll)-Pb 
interface. In order to describe properties related to the 
penetration of the wavefunctions into the substrate, the 
slab and infinite-square-potential-well models are inade- 
quate. 



3. Importance of quantum numbers 

The model based on the phase shifts or on the ef- 
fective width increase of the infinite potential well re- 
produces the eigenvalues and eigenfunctions of the self- 
consistent calculation. Because the self-consistently cal- 
culated eigenvalues agree with the experimental ones,^^ 
the latter can be fitted with 5 « 4ao. Nevertheless, Otero 
ef a/piS obtain a good agreement with 5 = 0. 

To explain the contradiction we have compared the 
infinite-square-potential-well energy spectrum (Eq. Q 
obtained with 5 = 0) and our model results correspond- 



ing to (5 ^ 0. In the latter the quantum number n for 
the eigenstates with the nearly constant eigenenergy of 
O.GbeV are n = 4, 7, 10, . . . for 2,4,6,. . . ML's of Pb, re- 
spectively. Actually, this is in agreement with the pseu- 
dopotential calculations for free-standing Pb slabs. But 
(5 = gives the corresponding quantum numbers as 
n = 3,6,9,... i.e., they are one unit smaller than the 
correct set. This explains why both models reproduce 
approximately the same experimental energy spectrum. 
The wavelength of the states n = 3, 6, 9, . . . at 0.65 eV 
in the infinite potential well is Ao.65 — 7.2ao. When 
we increase the width of the well by Ao.65/2 the states 
n = 4, 7, 10, . . . lie exactly at the same energy. This does 
not hold for energies far from 0.65eV and the eigenen- 
ergy spectra become remarkably different. Nevertheless, 
in the narrow energy window from -1 to 3 eV scanned in 
experiments, the energy spectra of both models are very 
similar and fit quite well the experimental results. As a 
matter of fact, an even better fit of experiments seems to 
be obtained by increasing the width of the potential well 
by Ao.65 a-iid using states with quantum numbers n two 
units larger than in the original infinite-square-potential- 
well model. 

In conclusion, the energy spectra can be fitted with 
several different thickness and quantum numbers. The 
danger is that other properties such as the wavefunctions 
or "magic" heights are wrongly determined. For exam- 
ple, the correct determination of the physical parameters 
reveals crucial in the description of the envelope function 
which plays a key role in determining the energetics of 
the quantum well states^ or the magnetic coupling in 
multilayer structures^ 



C. Determination of confinement barriers II 

1. Calculation of 5 from energy differences 

The disadvantage of the previous method for determin- 
ing the positions of the confining barriers is the necessity 
to measure a large number of QWS eigenenergies in or- 
der to label correctly the quantum states. However, there 
exists methods in which the knowlegde of the quantum 
number n is not needed MtSSiSLS^ We now elaborate this 
kind of methods. 

From the energy difference A„ between two consecu- 
tive states, one can obtain an effective thickness D' = 
D + 6a oi the infinite potential well as 



A.. 



7r2(2n-f 1) 



fcoTT 



(12) 



2{D + 6o)^ D + Sq 

Above, D is the thickness of the positive background den 
sity of the overlayer slab or, in the case of experiments 
the measured thickness, fco = 



_ 7r(rt+0.5) 
D+S„ 



is a wavevector 



depending on n and D'. As before, Sq is the difference 
between the measured nominal value and the effective 
one confining the electrons. The latter takes into ac- 
count, e.g., the effect of the electron spill out into the 



9 



80 



50 



30 



20 
















o 


uncorrected formula 


















o 


corrected formula 






+ 


exact 






























^ ' 














- 


















uncorrected 


8o=14.2ao 








corrected 


5o=5.40ao 








exact 


So =4.85 ao 



10 



20 



30 



40 



50 



60 



FIG. 6: Effective thickness fco7rA~^ calculated from the the- 
oretical QWS energies in Fig. (energy window from -1 to 
3.5 eV) as a function of the Pb jellium slab thickness. The 
diamonds and dashed line correspond to the data and fit ac- 
cording to Eq. 1)13^ . The circles and the solid straight line are 
the data and the fit according to Eq. 1151 1. The crosses are 
the exact effective thickness calculated from D' — 7m/\/2e„. 
Data for 23 ML's are used in the fitting, but only data for 11 
ML's are plotted for clarity. 
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FIG. 7: Energy eigenvalues (bottoms of subbands) of a free- 
standing 8 ML thick Pb slab as a function of the quantum 
number n. The open circles and the filled diamonds give 
the results of the stabilized-jellium and infinite-potential- well 
models, respectively. In the latter, the width of the well cor- 
responds to 8 ML widened by 4.5 ao in order to take into 
account the electron spill out. The vacuum level defines the 
zero energy level. The inset shows the energy differences be- 
tween the consecutive states as a function of the energy. The 
energy range of the QWS's in the Pb/Cu(lll) system is given. 



vacuum or into the substrate, the unknown thickness of 
the wetting layer, and effects due to stress or relaxation 
at the boundaries. 

We consider QWS's in overlayers of different thickness 
and within a given window, e.g. the QWS's shown in 
Fig. 13 We write Eq. (O as 



Then, by plotting A""'^ as a function of D and by fitting 
a straight line to the data, the slope gives fco averaged 
over the energy window in question. For example, the 
theoretical data of Fig. Ogives ko — O.SSqq^ which cor- 
responds to the kinetic energy of 10.5eT^, i.e. an energy 
'^l eV above the Fermi level. 

The original theoretical data and the fit are shown 
in Fig. by the diamond markers and by the dashed 
straight line, respectively. We plot D' — /cqI'A""'^, instead 
of A"-'^, as a function of D. This is because we want to 
include the exact values D' — nn/ ^/2en and the results 
of Eq. (|15|l (see below). The intersection of the dashed 
straight line with the vertical axis gives Sq = 14ao. This 
value is much larger than the value obtained for the total 
shift of the potential barriers in Section llll Bl Namely, at 
the energy of leV above the Fermi level we obtain from 
Fig. Ub) that Sq = 4.5ao. 

The method described above (Eq. (|13|l ) has been used 
to determine the number of Pb wetting layers in the STS. 
Our finding that the method overestimates the distance 
between the confinement barriers may explain the large 



^g^jyggiij20j2i^ obtained by this scheme in comparison 
with other experiments^*S*S The reason for this disagree- 
ment is the basic fact that Eqs. (|12|l and p3(l correspond 
to the infinite potential well, whereas a finite potential 
well is closer to reality. The infinite-square-potential- well 
model with the shifts 5q of the confining barriers repro- 
duces quite accurately the energy spectrum as is demon- 
strated in Fig. [3 The figure shows that the infinite 
potential well with the thickness of 8ML-|-4.5ao mimics 
the energy spectrum of a 8 ML unsupported Pb jellium 
slab. But the reasoning to the opposite direction, the 
determination of the Sq value from the real energy spec- 
trum cannot be done using the scheme. The explana- 
tion is that the energy eigenvalue differences A„ between 
the consecutive states behave differently in the infinite- 
square-potential- well and in real systems. The inset of 
Fig. 13 shows that the infinite potential well results in a 
monotonically increasing A„ whereas in the more real- 
istic stabilized-jellium model A„ decreases close to the 
vacuum level. The wrong trend in the infinite-square- 
potential-well model is compensated by the erroneously 
large effective width of the Pb slab or overlayer obtained 
by applying Eq. (fT^ . In the next subsections we suggest 
a method purified from this effect. 



2. Corrected formula for the calculation of So 

The option we choose to correct the overestimation in- 
herent in Eq. H13|) is to introduce the effect of the finite 
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potential barrier. This can be done by assuming that 
the surface shift is energy-dependent as (5 = + ^i^)- 
Here, 5q is the mean value we want to determine. Our 
aim is to obtain information about the electron confine- 
ment strength through the Sq parameter, which is energy- 
independent but it is reproduces satisfactorily the energy 
spectrum (see Fig. [7||. Nevertheless, it is necessary to use 
the energy-dependent 5{s) function to obtain relevant Sq 
values. 

The energy difference between the successive states can 
be obtained to the first order as A„ = 3^ increases by 
unity). The energy derivative for a given thickness D' is 



de TT^n ir'^n? 



(14) 



where 5{e) is the energy derivative of the surface shift and 
it depends on energy. This equation has to be evaluated 
for a given n. For the infinitely deep potential well <5 = 
and we do not recover the exact result of Eq. H12I) . This 
deficit can be corrected by evaluating the right-hand side 
of Eq. ifT^ at (n+1/2). The corresponding 5 is also eval- 
uated at (e„+i + £„)/2 using the self-consistent 5 values 
shown in Fig. Qs). Then, rearranging the terms of Eq. 
(|14|l we obtain the corrected formula for the thickness as 



TrfcoA: 



kl5 



D' = D + 5. 



(15) 



Here, it is accurate enough to use the feg obtained previ- 
ously with Eq. (|12|) . Omitting the dependence of 5 on 
energy on the right hand side of Eq. p5|l . i.e. 5 — S^^we 
can calculate the value of 6q by fitting a straight line for 
Z?' as a function of D. 

The circles and the solid straight line in Fig. |B| give 
the D' values and the fit, respectively, corresponding to 
Eq. H15|l . The new mean value = 5.4ao is consistent 
with the electron spill out shown in Fig. ^p) and it 
is much better than the value of 14ao obtained with the 
uncorrected formula (|13|l . Namely, the exact values D' ~ 
'nnj \J2en (shown as crosses in Fig. 0J| give = 4.85ao. 



3. Analytic models for the k^S correction 

In order to apply the corrected scheme of the previous 
subsection the energy-dependent derivative 6{e) has to 
be known. We study now the reliability of different ana- 
lytical models in its estimation. These models allow us to 
extend the previous analysis to other substrates or over- 
layers without doing self-consistent electronic structure 
calculations. 

The finite-square-potential-well model does not pro- 
vide a correction large enough, because the potential is 
not a continuous function of energy. The 6 values are too 
small (see Fig. Ela)). The image-potential model (Eq. 
(10)) and the empiric phase shifts (Eq. (|ll|l 'l provides 
results in a much better agreement with experiments. 
We want to emphazise that even if the analytical phase 
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FIG. 8: Derivative 5{e) of the effective thickness with re- 
spect to the energy according to different models. The solid 
and dash-dotted curves give the self-consistent results for the 
interface and surface components, respectively. The dotted 
curve shows the analytical result for the surface component 
obtained by using the image potential model of Eq. llOt . 
The dash-dotted curve gives the result for the interface com- 
ponent obtained by using the empirical formula (IIH . The 
energy origin is at the Fermi level. 



shifts may not produce the same Sq values as the self- 
consistent calculations (for example, due to the fact that 
the phase shift is determined with the module of 27r) they 
reproduce reasonably well the energy derivative 6 needed 
in the correction in Eq. H15(l . 

Fig. |S1 shows the derivatives S from the self-consistent 
calculation and from the analytical models of Sec. Ill CI 
The agreement is quite good, even if the derivatives of the 
self-consistent calculation are generally smaller. Compar- 
ing self-consistent and analythical curves we notice that 
small differences in the phases (Fig. (0))) produce big dif- 
ferences in Fig. |S1 As a matter of fact, at energies close 
to the vacuum the image-potential results could be more 
relevant than the self-consistent DFT results, which are 
affected by the exponential decay of the DFT potential. 
Because according to Fig. 8 the correction is large at 
energies close to the vacuum level and also close to the 
bottom of the projected energy gap of the substrate, we 
recommend the use of QWS's at intermediate energies 
when determining the overlayer thickness. 

Table HI shows a collection of Sq values obtained for 
several systems using different approximations. In addi- 
tion to our theoretical results for the Pb/Cu(lll) and 
the free-standing Pb slab systems we analyze our simi- 
lar results for Na/Cu(lll) and the experimental QWS 
spectrum of Pb/Cu(lll) [I^- The "exact <5o" values are 
obtained by fitting a straight line through the exact re- 
sults D' = Trn/ \J1en- In general, the uncorrected Eq. 
(|13|l gives for the Pb systems (5o values nearly three times 
larger than the exact one. In contrast, the corrected value 
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TABLE L So's (in Bohr radii ao) determining the effective widths of the infinite potential weUs fitting the QWS energy spectra. 
The values are calculated from different theoretical and experimental spectra by using the uncorrected sheme of Eq. 113^ as 
well as the scheme of Eq. 1151 ') with the numerical correction obtained from the self-consistently calculated results or from 
the analytical image-potential and empiric formulae. The exact values are obtained by employing the equation D' = Trn/y'2e„ 
with the known n quantum numbers. The maximun overlayer thickness (in ML's) used in fitting are given in the second 
column. The numbers in parenthesis give the linear regression errors. The QWS of the DFT calculations used in fitting span 
the energy range from ~ leV below Fermi level, i.e. from the bottom of the Cu(lll) band gap, to ~ leV below vacuum level 
(for Pb/Cu(lll) corresponds to the experimentaly measured energy range of QWS's). 



System 


Up to 


<5o from Eq. itHl 
uncorrected 


So from Eq. IT^ 
numerical correction 


So from Eq. itTCl 
analytical correction 


Exact So 


Free-standing Pb 


10 ML 


12(3) 


5.5(0.6) 


6(3) 


3.9(0.3) 


Pb/Cu(lll) 


23ML 


12(1) 


5(2) 


3(2) 


4.7(0.1) 


Pb/Cu(lll) 


lOML 


13(1) 


6(1) 


4(2) 


4.8(0.3) 


Experiments— 


24ML 


15(5) 


7(5) 


5(5) 


4.0(0.2) 


Na/Cu(lll) 


20ML 


7(3) 


1(3) 


3(4) 


1.2(0.3) 



offers a sub-monolayer accuracy. 

For the Na/Cu(lll) system the difference between the 
corrected and uncorrected <5o is large as well. Although 
the wavelength of Na is larger than for Pb, Sq does not 
scale with it. The deeper potential of Cu(lll) than 
the Na potential (the opposite to the Pb/Cu(lll) sys- 
tem) shifts the barrier inside the Na, i.e. the value of 
'^Cti(iii)-JVa is negative. This explains the small value of 
So- 

It is noticeable that applying Eq. 1)15(1 with the nu- 
merical and the analytical corrections to the experimen- 
tal results of ref. [l^ we obtain Sq — 7.3ao and 4.7ao, 
respectively. These values are much smaller than the 
uncorrected value of 14.9ao. When we compare with the 
exact value of 4.0ao, which is obtained by identifying first 
the n quantum numbers. The linear regression error is 
anyway substantial. In addition to the errors due to our 
first-order approximation in Eq. (14), the experimental 
dispersion in the QWS spectrum cause uncertainty. In 
general, for self-consistent calculations the numerical cor- 
rection provides better determination of Sq compared to 
the analytical one. But in the case of the experimental 
data the analytical correction is better. 

When analyzing experiments, the So value obtained by 
the corrected scheme of Eq. (15) can be used as an initial 
parameter to determine the n quantum number. Then, 
employing the exact equation to fit the energy spectrum, 
improved results are obtained. 



IV. CONCLUSIONS 

We have performed self-consistent DFT calculations to 
study the confinement barriers of electrons in Pb islands 
grown on the Cu(lll) substrate. Calculations have been 
done for free-standing Pb slabs and Pb slabs on Cu(lll). 
Pb was described by stabilized jellium and the Cu(lll) 
substrate by a ID-pseudopotential. The model repro- 
duces the most important physical properties and gives 



results in a good agreement with experiments. 

The energies and wavefunctions of the quantum well 
states in the Pb slabs characterize the confinement bar- 
riers at the Pb- vacuum surface and at the Cu(lll)-Pb in- 
terface. We have analyzed these states by using the phase 
accumulation model and by determining effective widths 
of infinite potential wells reproducing the energies. The 
Pb- vacuum phase shift is in a good agreement with more 
realistic pseudopotential calculations. The Cu(lll)-Pb 
phase shift or the effective width of the potential well ac- 
counts for the confining strength of the Cu(lll) energy 
gap. This strength is weaker than that of the Pb- vacuum 
barrier. 

The information provided by our calculations and anal- 
ysis allows to improve the interpretation of QWS spec- 
tra measured by scanning tunneling spectroscopy. More 
specifically, we have shown that a formula commonly 
used in the literature results easily in the overestima- 
tion of the effective width of the infinite potential well. 
We have offered an alternative expression to correct that 
deficiency, specially important for electronic high density 
(small Xf) metals. The obtained results can be used to 
estimate the width of the potential well and to determine 
the quantum numbers for a more accurate analysis of the 
confinement barriers. 

Finally, the ID-pseudopotential scheme provides a 
method for future studies of nanostructures on solid sur- 
faces. Despite this work is focused on STS experiments 
and Pb/Cu(lll) system, our results can be qualitatively 
applied to metallic nanoislands on semiconducting sub- 
strates, e.g. Pb or Ag on Si, as well as to photoemission 
experiments. 
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